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:【 目 的 】] 肚 贝 科 食性 高 度 分 化 。 本 研究 旨 在 通过 建立 七 星球 虫 Coccinella septempunctata 触角 
触角 嗅觉 相关 基因 与 食性 分 化 的 关系 。【 方 法 】 采 用 Ilumina HiSeq 4000 高 
人 肚 虫 成 下 触角 转录 组 进行 测序 、 组 装 、 注 释 , 挖 握 嗅 觉 相关 基因 ,并 与 已 发 表 的 
tA ZJ k Henosepilachna viginyioctopunctata 转录 组 进行 比较 。 【结果 】 共 获得 七 星球 xh f 
转录 组 31 775 条 unigenes, 其 中 69.71% 的 序列 得 以 注释 ,NR 数据 库 中 注释 最 多 ,为 20 539 条 。 
注释 信息 ,挖掘 到 27 个 嗅觉 相关 基因 ,包括 1 个 气味 结合 蛋白 (odorant binding protein, OBP) 基 
13 个 化 学 感受 蛋白 (chemosensory protein, CSP) 基 因 ,4 个 气味 受 体 (odorant receptor, OR ) 基 
个 味觉 受 体 (gustatory receptor, GR) 基因 和 2 个 感觉 神经 元 膜 蛋 和 白 (sensory neuron vnd 
protein, SNMP) X, 483] E35 , AE RE 35 — FE SIC A Asp XE 8138 个 嗅觉 相关 基因 , 包 
4&-5 X SK 束 中 未 发 现 的 1 工 个 离子 型 受 体 (ionotropic receptor，IR) 基 因 。 在 各 类 型 嗅觉 相关 基因 中 ， 
Ao — TAE SK de RKA OBP 基因 比例 (13.16% ) 高 于 七 星 靳 虫 触角 转录 组 的 (3.70% ) ,而 七 星 
上 黑 忠 触角 转录 组 的 GR 基因 比例 (25.93% ) 则 高 于 茄 二 十 八 星 靳 虫 转录 组 的 (13.16% )。[【 结论 】 
触角 嗅觉 相关 基因 数目 不 是 昆 求 食性 分 化 的 主因 。 REIS eu 触角 转录 组 学 资源 , 初 
步 探 讨 了 嗅觉 相关 基因 同 默 虫 食 性 分 化 的 关系 ,为 了 解 球 虫 帮 至 昆虫 食性 分 化 的 分 子 基 础 提供 了 
信息 。 
关键 词 : GEMER; 触角 转录 组 ; 食性 分 化 ; 基因 注释 ; 嗅觉 相关 基因 
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Abstract: [ Aim] Coccinellidae show highly diverse feeding. This study aims to explore the relationships 
between olfaction-related genes in antennae and feeding divergence based on the establishment of the 


antennal transcriptome of Coccinella septempunctata. [Methods] The antennal transcriptome of C. 
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septempunctata adults was sequenced using an Illumina HiSeq 4000 platform and subjected to assembly, 
annotation and comparative analysis with the previously published transcriptomes of Henosepilachna 
viginyloctopunctata , and the olfaction-related genes were identified. [Results] A total of 31 775 unigenes 
were obtained in C. septempunctata antennal transcriptome, 69. 7196 of which were annotated, and the 
highest number of unigenes (20 539) were annotated in the NR database. According to the annotation 
information, 27 olfaction-related genes including 1 odor binding protein ( OBP) gene, 13 chemoreceptor 
protein ( CSP) genes, 4 odorant receptor ( OR) genes, 7 gustatory receptor ( GR) genes and 2 sensory 
neuron membrane protein ( SNMP) genes were identified. In parallel, 38 olfaction-related genes were 
identified in the transcriptomes of H. viginyioctopunctata , including 1 ionotropic receptor (IR) gene that 
was not found in C. septempunctata. Among the various types of olfaction-related genes, the proportion 
(13. 1696 ) of OBP genes in the transcriptomes of H. viginyioctopunctata was higher than that (3. 7096 ) 
in the antennal transcriptome of C. septempunctata, while the proportion (25. 9396 ) of GR genes in the 
antennal transcriptome of C. septempunctata was higher than that (13. 1696 ) in the transcriptomes of H. 
viginyioctopunctata. [ Conclusion] The number of olfaction-related genes in the antennae may not make a 
predominant contribution to insect feeding divergence. This study obtained the antennal transcriptomic 
data of C. septempunctata and preliminarily explored the relationships between olfaction-related genes and 
feeding divergence in ladybirds, providing information for understanding the molecular basis of feeding 
divergence in ladybirds and even in insects. 


Key words: Coccinella septempunctata; antennal transcriptome; feeding divergence; gene annotation; 


olfaction-related gene 








嗅觉 是 昆虫 食物 选择 求偶、 交配 及 产 卵 等 各 种 
行为 的 基础 ,对 其 种 群生 存 、 繁 衍 有 重要 作用 (Field 
et al., 2001; Sato and Touhara, 2009; Zhang et al., 
2015; Gadenne et al., 2016)。 昆 虫 触角 上 着 生 许 多 
嗅觉 感 器 上 且 有 各 种 嗅觉 蛋 白 ,不同 的 感 器 和 嗅觉 蛋 
白 可 识别 不 同 的 气味 物质 进而 调控 昆虫 的 行为 活动 
( Brown et al., 2013; 郑海霞 等 , 2018 ) 。 目 前 ,昆虫 
触角 中 研究 最 为 广泛 的 嗅觉 相关 和 蛋白 主要 包括 : 气 
味 结合 蛋白 (odorant binding protein, OBP) 、 化 学 感 
ZE H (chemosensory protein, CSP) 、 嗅 党 /气味 受 
体 ( olfactory odorant receptor，OR ) 、 离 子 型 受 体 
(ionotropie receptor, IR )、 味 觉 受 体 (gustatory 
receptor, GR ) 和 昆虫 感觉 神经 元 膜 蛋白 (sensory 
neuron membrane protein, SNMP) 六 大 类 ( Leal, 
2013; Zhang et al., 2016; Wang J et al., 2017) , ?& 
Bii , Ak F Eb B f £e Pe e 2H ^E TZ i ES TH OR AE DAT D. 28 
被 广泛 开展 。Zhang 5$ (2016) A 3€ E AIR 
Hyphantria cunea 触角 转录 组 中 鉴定 了 124 个 嗅觉 
相关 基因 ,主要 包括 OR, OBP 和 CSP 基因 等 ,并 分 
析 了 OBP 和 CSP 基因 的 表达 模式 。Jia 等 (2016) 采 
用 了 相似 的 研究 策略 鉴定 了 桃 星 量 Conogethes 
punctiferalis 瞧 雄 成 虫 触角 中 的 嗅觉 相关 基因 ,并 对 
呈 性 别 差 异 表达 的 嗅觉 相关 基因 进行 了 鉴定 。 有 研 























窒 从 云南 同 域 生境 下 云 商 松 毛 虫 Dendrolimus houi 
BIMBE D. kikuchü 触角 转录 组 中 鉴定 了 一 批 
嗅觉 相关 基因 ,进行 了 属 内 和 目 间 比较 和 基因 表达 
水 平分 析 , 探 讨 了 基因 功能 的 潜在 相似 性 (Zhang et 
al., 2014)。 然 而 ,迄今 为 止 未 见 从 嗅觉 相关 基因 的 
角度 探讨 其 在 不 同 食性 昆虫 中 的 差异 。 昆 虫 触角 嗅 
觉 相 关 基 因 是 否 是 其 食性 分 化 的 主要 因素 ,昆虫 嗅 
觉 相关 基因 数量 及 不 同 食性 昆虫 间 基因 数量 差异 对 
了 解 昆 虫 间 食性 分 化 能 提供 什么 信息 等 ,都 是 有 待 
探讨 的 问题 。 
SL di Æ H xH H ( Coleoptera ) JE R F} 
( Coccinellidae ) 昆虫 的 总 称 ,涵盖 90 个 属 超过 1 000 
个 种 ( Escalona et al., 2017), E: FP E E 3 E 
sk m OS dm ow B 
( Coccinellinae) ,是 以 蚜虫 为 主食 的 捕食 性 昆虫 
(Van der Werf et al., 2000) 。 飘 虫 科 昆虫 已 经 被 广 
泛 地 用 于 科学 研究 ,主要 集中 在 食性 可 塑性 的 发 育 、 
系统 发 育 .种 群 遗 传 、 峭 记 颜 色 多 样 化 和 捕食 性 (天 
政 ) 球 虫 对 其 他 害虫 的 生物 防治 等 (Escalona et al., 
2017; Matsubayashi et al., 2017; Saito et al., 2017; 
Tang et al., 2017) 。 值 得 注意 的 是 , 科 内 高 度 的 食性 
分 化 是 球 虫 科 的 一 大 特征 ( Escalona et al., 2017) ， 
大 部 分 标 虫 都 是 以 捕食 蚜虫 ( 如 棉 是 Aphis gossypii, 



































Coccinella septempunctata 





























6 期 杨 雪 娇 等 : EE B fü ff Pe oe 2H LS DEUS RH OS EDS 3p 719 











大 豆 蚜 A. glycines, $ BB Lipaphis erysimi 5g) W 
TER BPA H 4] rl, IAKA RETER E 
又 虫 。 男 外 ,还 有 少 部 分 种 类 为 植 食性 又 虫 , 主 要 包 
括 草食 性 的 食 植 又 虫 亚 科 ( Epilachninae) ,以 花粉 和 
蔬菜 等 为 食 ,多 数 为 农业 害虫 ,如 茄 二 十 八 星球 虫 
2002 ; 





Henosepilachna viginyioctopunctata 等 ( 庞 虹 ， 
Zhang et al., 2018) 。 

为 了 探索 昆虫 触角 嗅觉 相关 基因 在 昆虫 食性 分 
化 中 的 地 位 及 角色 ,本 研究 以 捕食 性 七 星球 虫 为 研 
究 对 象 , 运 用 Ilumina HiSeq 高 通 量 测序 技术 对 其 触 
角 转 录 组 进行 测序 ,采用 生物 信息 学 方法 对 其 基 
功能 注释 ,并 联合 我 们 以 前 测序 的 植 食性 茄 二 十 八 
星 昧 虫 转录 组 鉴定 并 比较 各 类 型 嗅觉 相关 基因 占 咱 
觉 相关 基因 总 数 的 比例 。 研 究 结果 将 加 深 对 球 虫 食 
性 分 化 分 子 基 础 的 了 解 , 也 为 驹 虫 科 天 敌 利 用 和 害 
虫 防治 提供 理论 依据 。 





























1 材料 与 方法 


1.1 WH 

E EBUR MEER RR AS PEJE B T8 KALET 
(海拔 10 m, 109°08'E, 21?25'N) , 。 为 了 避免 局 部 
环境 的 差异 ,七 星际 虫 分 别 采集 于 不 同 的 3 个 点 ， 
个 点 相距 至 少 30 km 以 上 。 因 七 星球 虫 成 虫 特征 明 
T ,种 类 鉴定 简单 易 行 ,野外 鉴定 带 回 实验 室 后 , 选 
取 活 动能 力 较 强 的 七 星际 虫 成 虫 200 头 , 切 下 触角 ， 
经 液 所 速冻 ,后 保存 于 - 807€ 超低温 冰箱 (Thermo 
Scientific, 美国 ) 待 用 。 
1.2 总 RNA 提取 

采用 RNeasy Mini Kit( 天 根 , 北京 ) 试 剂 盒 分 别 
提取 供 试 肉 雄 七 星球 虫 成 虫 各 100 头 触 角 总 RNA, 
经 NanoDrop@ ND1000 ( Thermo Scientific， 美 国 ) 核 
酸 浓度 测定 仪 检测 RNA 的 浓度 和 纯度 ,并 采用 
1.596 的 琼脂 糖 凝 胶 电 泳 和 Agilent 2100 ( 安捷伦 科 
技 , 美国 ) 精 确 检测 RNA 的 完整 性 。 等 量 混合 合格 
HWER- EDULE AfA A RNA 为 一 个 混 池 , 用 于 后 续 
转录 组 测序 。 
1.3 cDNA 文库 的 构建 和 测序 

合格 的 七 星 蒜 虫 触 角 总 RNA 使 用 DNase I 
(Qiagen， 德 国 ) 消化 残余 基因 组 DNA 后 ,按照 
Preparing Samples for Sequencing of mRNA 试剂 盒 说 
WE (lumina, 美国 ) 构 建 测序 文库 。 具 体 步 又: 用 
带 有 Oligo( dT) 的 磁 珠 寅 集 总 RNA 池 中 的 mRNA ; 
加 入 打上 断 试 剂 后 于 Thermomixer C ( Eppendorf, 4E 






































国 ) 中 适 温 将 mRNA 打 断 成 短片 段 (200 ~400 bp), 
以 打 断 后 的 mRNA 短片 段 为 模板 合成 cDNA 第 1 
链 , 随 后 在 DNA polymerase I 体系 下 ,以 dNTPs 为 原 
料 合成 cDNA 第 2 条 链 ,并 使 用 试剂 盒 纯 化 回收 、 粘 
性 未 端 修复 .cDNA 的 3 末端 加 上 碱 基 A 并 连接 接 
头 , 然 后 进行 片段 大 小 选择 ,最 后 进行 PCR 扩 增 ; 构 
建 好 的 文库 用 Agilent 2100 Bioanalyzer (安捷伦 科 
技 , 美国 ) 和 ABI StepOnePlus Real-Time PCR System 
( Applied Biosystems ,美国 ) 质 检 , 合 格 后 使 用 HiSeq 
4000 高 通 量 测序 平台 (Ilumina,， 美国 ) 的 PE150 bp 
双向 测序 技术 进行 测序 。 
1.4 转录 组 测序 数据 组 装 及 完整 性 评估 

利用 Trimmomatic vO. 39 ( Bolger et al., 2014) 软 
件 对 测序 得 到 的 原始 数据 (raw data) 进行 过 滤 ,去 
除 接头 序列 以 及 低 质量 reads ,获得 高 质量 的 数据 
(clean data) ,具体 步骤 如 下 : (1) 去除 包含 接头 的 
reads( 接 头 污染 );(2) 去 除 未 知 碱 基 N 含量 大 于 
5% 的 reads;(3 ) 去除 低 质量 的 reads ( 质量 值 低 于 15 
的 碱 基 占 该 reads 总 碱 基数 的 比例 大 于 30% 的 reads 
为 低 质 量 的 reads ) 。 利 用 Trinity 软件 对 clean data 
进行 序列 拼接 与 组 装 。 然 后 ,去 除 长 度 小 于 150 bp 
的 序列 ,并 通过 CD-HIT 软件 包 (http: // www. 
bioinformatics. org/ cd-hit/ ) 中 的 CD-HIT-EST 子 程序 
对 unigenes 以 95% 的 相似 性 进行 聚 类 ,并 选取 最 长 
的 序列 为 每 一 个 聚 类 艇 的 代表 ,以 去 除 unigenes 的 
宛 余 ,最 终 获 得 高 质量 的 七 星 标 虫 触角 unigenes 集 。 
使 用 BUSCO ( Benchmarking Universal Single Copy 
Orthologs) v3 软件 于 默认 参数 下 ,将 最 终 组 装 获得 的 
unigenes 通过 与 insecta odbO 数据 库 (https: // 
busco. ezlab. org/) 1 658 个 昆虫 进化 保守 的 直系 同 
源 基因 进行 比较 ,检测 每 个 unigene 的 完整 性 ,以 评 
估 转 录 组 组 装 质量 。 
1.5 转录 组 数据 分 析 及 嗅觉 相关 基因 鉴定 

将 1.4 节 所 得 unigenes, 使 用 本 地 Blastx 软件 分 
35 NR, NT, COG, Swiss-Prot 和 KEGG 等 数据 库 
进行 比 对 (FE 值 <1 x 10 7) ,获得 相应 的 注释 信息 。 
根据 NR 注释 信息 获得 基因 功能 信息 ,然后 提取 每 
个 基因 相应 的 GO 注释 信息 ,并 基于 GO 注释 把 基 
因 划 分 到 生物 学 过 程 、 细 胞 组 分 和 分 子 功能 等 三 亚 
类 。 利 用 COG 数据 库 注释 信息 对 直系 同 源 蛋 白 家 
族 进 行 功能 分 类 。KEGG 注释 信息 用 以 呈现 基因 参 
与 的 代谢 通路 。 茄 二 十 八 星球 虫 高 质量 的 unigenes 
集 及 相关 基因 功能 注释 信息 从 我 们 以 前 的 数据 集 相 
关 研 究 中 获取 (Zhang et al., 2018) 。 最 后 ,采用 以 
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往 文献 使 用 的 方法 (Jia et al., 2016; Zhang et al., 
2016; 郑海霞 等 , 2018 ) , 即 根据 基因 注释 信息 鉴定 
嗅觉 相关 基因 。 
2 结果 
2.1 七 星球 虫 触角 转录 组 序列 分 析 与 组 装 

七 星 标 虫 触角 转录 组 测序 获得 10.73 Gb 原始 
数据 。 质 控 后 获得 70 291 132 条 共 10. 54 Gb 高 质 
量 clean reads 数据 ,Q20 和 Q30 值 碱 基 比 例 分 别 为 
97.57% 和 93. 69% 。 利 用 Trinity 软件 对 获得 的 高 
质量 reads 进行 组 装 , 共 获得 46 687 条 contigs ,总 长 
度 为 41 887 759 bp ,平均 长 度 为 897 bp,N50 的 长 度 
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为 1 871 bp, N7O 的 长 度 为 1 061 bp, N90 的 长 度 为 
311 bpo Xf contigs 进一步 进行 聚 类 和 去 宛 余 ,得 到 
31 775 条 unigenes ,总 长 度 为 34 872 097 bp ,平均 长 
1 097 bp, N50 的 长 度 为 2 010 bp; GC 含量 为 
38. 99% , unigenes 长 度 分 布 (图 1) 显示 ,在 长 度 
1 kb 以 上 的 unigenes 为 13 107 条 , 占 unigenes 总 数 
的 41. 25%。BUSCO 分 析 结 果 显 示 83. 22% 的 
unigenes 为 完整 基因 ,7. 13% 是 部 分 片段 , 剩 下 的 
unigenes 没有 比 对 到 任何 序列 。 为 了 使 相关 学 者 使 
用 方便 ,避免 前 期 数据 组 装 的 不 便 ,unigenes 集 已 释 
放 到 公共 存储 数据 库 Figshare ( https: // doi. org/ 
10. 6084/m9. figshare. 10259660. v1) 供 直 接 下 载 
使 用 。 
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图 1 七 星球 虫 触角 转录 组 unigenes 长 度 分 布 











Fig. 1 


2.2 七 星球 虫 触角 转录 组 基因 功能 注释 

七 星 标 虫 触角 转录 组 中 ,共有 22 149 条 
(69.71% ) unigenes 获得 注释 。 其 中 NR 数据 库 注 
释 到 的 unigenes 数目 最 多 ,为 20 539 条 (64. 6496 ) ; 
GO 数据 库 注释 的 unigenes 最 少 , 为 4 863 条 
(15. 30% ) ;NT 数据 库 注释 到 7 576 2& (23. 8496) ; 
Swiss-Prot 数据 库 注 释 到 16 862 条 (53. 0796) ; 
KEGG 数据 库 注释 到 16 445 2& (51. 7596 ) ; COG 数 
据 库 注释 到 9 342 2& (29. 4096 ) 。 

基于 NR 注释 信息 ,七 星 标 虫 触角 转录 组 中 共 
注释 到 20 539 条 unigenes, 其 中 与 赤 拟 谷 资 
Tribolium castaneum 的 相似 基因 最 多 ,所 占 比 例 为 
62.70% ,其 后 依次 为 中 欧 山 松 大 小 震 Dendroctonus 
ponderosae (7. 79% )、 丽 量 师 集 金 小 蜂 Nasonia 





Length distribution of all assembled unigenes in the antennal transcriptome of Coccinella septempunctata 


vitripennis ( 1. 2196 ) , Sij X3 WW Acyrthosiphon pisum 
(1.05% ) ,其 他 物种 占 27.26% (图 2) 。 基 于 GO 数 
据 库 对 七 星球 虫 触 角 转 录 组 进行 功能 分 类 统计 (图 
3 ) , 共 得 到 22 020 条 功能 注释 信息 ,分布 于 4 863 条 
unigenes ,生物 学 过 程 (9 841 条 unigenes) 中 代谢 过 
程 基因 最 多 (2 373 条 unigenes) ,其 次 是 细胞 进程 
(2 325 条 unigenes ) 和 单 生 物 进 程 (1 695 条 
unigenes) , 其余 的 均 低 于 1 000 条 unigenes ; 细胞 组 
分 (6 406 条 unigenes) 中 细胞 部 分 和 细胞 最 多 , 均 为 
1 454 条 unigenes ,其 次 是 细胞 器 (931 条 unigenes) 、 
膜 ( 800 条 unigenes )、 大 分 子 复合 物 (641 条 
unigenes)、 膜 部 分 (529 条 unigenes ) , 其 余 均 低 于 
500 条 unigenes; 分 子 功能 (5 773 条 unigenes ) 中 催化 
活性 有 2 430 条 unigenes, 结合 有 2 350 条 unigenes, 
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其 余 均 低 于 500 条 unigenes。 
27.26% 在 COG 数据 库 中 ,9 342 条 unigenes 被 分 为 25 
个 大 类 , 共 包 含 14 906 条 注释 信息 (图 4)。 其 中 数 
量 最 多 的 是 仅 一 般 功 能 预测 类 (2 841 条 ， 
19. 06% ) ,其 余 依 次 为 翻译 、 核 糖 体 结构 和 生物 合 
1.05% 成 类 (1 197 %, 8.0396) , 复制 、 重组 和 修复 类 
人 (1 136 条 , 7.62% ) 等 。 根 据 KEGG 数据 库 对 七 星 
7395 ni 球 虫 触角 转录 组 的 注释 信息 ,16 445 条 unigenes 注 
释 到 31 645 个 代谢 通路 ( 表 1) ,有 9 个 代谢 通路 包 
含 300 条 以 上 的 unigenes: 代 谢 通路 (649, 3.95% ) 、 
LOO AR SPA UA D ELaR M ee 
IRUIN Nasonia vitripennis 2.39% )、 艾 巴 氏 病毒 感染 (382,， 2.32% ) RISA 
Vor A M (370, 2.25% ) MEUIE (368, 2.2490) ,癌症 中 的 蛋 
图 2 七 星 蒜 下 触角 unigenes 在 NR 数据 库 中 的 物种 分 布 。 ”日 多 糖 (362, 2. 20% ) RNA 转运 (328, 1.99%), 














Fig. 2 Species distribution of antennal unigenes of 
Coccinella septempunctata in the NR database 


行为 Behavior 

生物 粘 附 Biological adhesion 

生物 相 Biological phase 

生物 调控 Biological regulation 

细胞 杀伤 Cell killing 

细胞 组 分 组 成 或 生物 合成 Cellular component organization or biogenesis 
细胞 进程 Cellular process 

发 育 进 程 Developmental process 

生长 Growth 

免疫 系统 进程 Immune system process 

定位 Localization 

移动 Locomotion 

代谢 过 程 Metabolic process 

多 细胞 生物 进程 Multicellular organismal process 

多 生物 过 程 Multirorganism process 

生物 过 程 负 调控 Negative regulation of biological process 

生物 过 程 正 调控 Positive regulation of biological process 

生物 过 程 调控 Regulation of biological process 

生殖 Reproduction 

生殖 过 程 Reproductive process 

应 激 响 应 Response to stimulus 

信号 Signaling 

单 生物 进程 Single-organism process 

细胞 Cell 

细胞 连接 Cell junction 

细胞 部 分 Cell part 

胶原 三 聚 体 Collagen trimer 

胞 外 基质 Extracellular matrix 

胞 外 区 域 Extracellular region 

胞 外 区 域 部 分 Extracellular region part 

大 分 子 复合 物 Macromolecular complex 

膜 Membrane 
膜 封闭 腔 Membrane-enclosed lumen 

膜 部 分 Membrane part 

拟 核 Nucleoid 

细胞 器 Organelle 

细胞 器 部 分 Organelle part 

突 触 Synapse 

突 触 部 分 Synapse part 

病毒 粒子 Vrion 

病毒 粒子 部 分 Virion part 

抗 氧化 活性 Antioxidant activity 

结合 Binding 

催化 活性 Catalytic activity 

电子 载体 活性 Electron carrier activity 

酶 调控 活性 Enzyme regulator activity 

鸟 彰 酸 交换 因子 活性 Guanyl-nucleotide exchange factor activity 
金属 伴侣 活性 Metallochaperone activity 

分 子 转 导 器 活性 Molecular transducer activity 

结合 转录 因子 活性 Nucleic acid binding transcription factor activity 
转录 因子 活性 Protein binding transcription factor activity 
蛋白 标签 Protein tag 

受 体 活性 Receptor activity 

结构 分 子 活性 Structural molecule activity 
翻译 调控 器 活性 Translation regulator activity 
转运 体 活 性 Transporter activity 















































肌 动 蛋白 骨架 调节 (322, 1.96% ) ;包含 150 ~ 300 
条 unigenes 的 代谢 通路 有 18 个 ;其 余 代 谢 通路 所 含 
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Fig. 3 GO function annotation of antennal unigenes of Coccinella septempunctata 
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翻译 、 核 糖 体 结构 和 生物 合成 Tanslation, ribosomal structure and biogenesis 
转录 Transcription 
信号 转 导 机 制 Signal transduction mechanisms 
次 生 代谢 物 生物 合成 、 转 运 和 催化 Secondary metabolites biosynthesis, transport and catabolism 
RNA 进 程 和 修饰 RNA processing and modification 
复制 、 重 组 和 修复 Replication, recombination and repair 
翻译 后 修饰 、 蛋 白 转换 、 伴 但 Posttranslational modification, protein turnover, chaperones 
核 音 酸 转运 与 代谢 Nucleotide transport and metabolism 
核 结构 Nuclear structure 
转运 和 代谢 Lipid transport and metabolism 
细胞 内 运输 、 分 泌 和 宫 泡 转运 Intracellular trafficing, secretion, and vesicular transport 
无 机 离子 转运 与 代谢 Inorganic ion transport and metabolism 
仅 一 般 功 能 预测 General function prediction only 
未 知 功能 Function unknown 
细胞 外 结构 Extracellular structures 
能 量 生 产 和 转换 Energy production and conversion 
防御 机 制 Defense mechanisms 
细胞 骨架 Cytoskeleton 
辅酶 转运 与 代谢 Coenzyme transport and metabolism 
染色 质 结构 与 动力 学 Chromatin structure and dynamics 
细胞 壁 / 膜 / 包 膜 生物 发 生 Cell wall/membrane/envelope biogenesis 
细胞 死亡 Cell motility 
细胞 周期 控制 、 细 胞 分 裂 、 染 色 体 分 型 Cell cycle control, cell division, chromosome partitioning 
碳水 化 合 物 运输 和 代谢 Carbohydrate transport and metabolism 
氨基 酸 转运 和 代谢 Amino acid transport and metabolism 
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Fig. 4 COG function classification of antennal unigenes of Coccinella septempunctata 


R1 七 星球 虫 触角 转录 组 KEGG 通路 分 析 


Table 1 KEGG pathways in the antennal transcriptome of Coccinella septempunctata 
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通路 ID 代谢 通路 基因 数量 ii EE C96 ) 
Pathway ID Metabolic pathways Number of genes Proportion 
ko01100 代谢 通路 Metabolie pathways 649 3.95 
ko05200 癌症 中 通路 Pathways in cancer 538 3.27 
ko04144 WAE Endocytosis 393 2.39 
ko05169 艾 巴 氏 病毒 感染 Epstein-Barr virus infection 382 2.32 
ko00230 喇叭 代谢 Purine metabolism 370 2.25 
ko04145 吞噬 体 Phagosome 368 2.24 
ko05205 癌症 中 蛋白 多 糖 Proteoglycans in cancer 362 2.20 
ko03013 RNA 转运 RNA transport 328 1.99 
ko04810 肌 动 蛋白 骨架 调节 Regulation of actin cytoskeleton 322 1.96 
ko04151 PI3K-Akt 信号 通路 PIBK-Akt signaling pathway 265 1.61 
ko03010 核糖 体 Ribosome 264 1.61 
ko04510 粘 附 斑 Focal adhesion 262 1.59 
ko01200 碳 代谢 Carbon metabolism 258 1.57 
ko03040 剪接 体 Spliceosome 253 1.54 
ko04015 Rapl 信号 通路 Rapl signaling pathway 251 1.53 
ko05016 享 廷 顿 氏 病 Huntington’ s disease 248 1.51 
ko04024 cAMP 信号 通路 cAMP signaling pathway 207 1.26 
ko04919 甲状 腺 激素 信号 通路 Thyroid. hormone signaling pathway 286 1.74 
ko05130 致 病 性 大 肠 杆菌 感染 Pathogenic Escherichia coli infection 282 1.71 
ko04360 轴 突 导向 Axon guidance 281 1.71 
ko03015 mRNA 监控 通路 mRNA surveillance pathway 276 1.68 
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续 表 1 Table 1 continued 





























通路 ID 代谢 通路 基因 数量 ii EE C6 ) 
Pathway ID Metabolic pathways Number of genes Proportion 
ko05100 表皮 细胞 的 细菌 入 侵 Bacterial invasion of epithelial cells 275 1.67 
ko05166 HTLV-I 感 当 HTLV-I infection 174 1.06 
ko04141 内 质 网 蛋白 质 加 工 Protein processing in endoplasmic reticulum 172 1.05 
ko04611 血小板 活化 Platelet activation 169 1.03 
ko05131 志 贺 氏 菌 病 Shigellosis 165 1.00 
ko05010 阿尔 茨 海 默 氏 病 Alzheimer’ s disease 163 0.99 
其 他 通路 Other pathways <150 














unigenes 数量 均 低 于 150 条 。 
2.3 七 星球 虫 触角 转录 组 嗅觉 相关 基因 分 析 
基于 六 大 公共 数据 库 得 到 的 unigenes 注释 信 
息 ,27 个 噢 觉 相 关 基 因 在 七 星 飘 虫 触角 转录 组 中 被 
鉴定 ( 表 2): 1 个 气味 结合 蛋白 (OBP) 基因 
(3.7096), 13 个 化 学 感受 蛋白 (CSP) 基因 
(48. 15% ) .4 个 气味 受 体 (odorant receptor, OR) 基 
因 (14. 48% ) .7 个 味觉 受 体 (GR ) 基因 (25. 93% ) 
和 2 个 感觉 神经 元 膜 蛋白 (SNMP) 基因 (7.41% ) 。 


另外 ,在 茄 二 十 八 星 球 虫 转录 组 中 鉴定 到 38 个 嗅觉 
相关 基因 ,包括 5 个 OBP 基因 (13.16% ) , 17 个 
CSP 基因 (44.74% ), 4 个 OR 基因 (10.53%),1 
A IR 基因 (2.63% ), 5 个 GR 基因 (13.16% ) 和 6 
个 SNMP 基因 (15. 7996) 。 比 较 发 现 , 茄 二 十 八 星 
昧 虫 转录 组 各 类 型 嗅觉 相关 基因 中 ,OBP 基因 比例 
(13.16% ) 高 于 七 星球 虫 触角 转录 组 的 (3.70% ) ,而 
七 星球 虫 触角 转录 组 中 GR 基因 的 比例 (25. 93% ) 高 
Toti — T /NESIUR (I3. 1696) 。 


表 2 LESE H ERARA BAIE 


Table 2 Olfaction-related genes in the transcriptomes of Coccinella septempunctata and 


Henosepilachna viginyioctopunctata 



















































































EC B fh f FEH Jii — FN EE BUB SEE 
基因 Antennal transcriptome of C. septempunctata Transcriptomes of H. viginyioctopunctata 
基 
Genes 嗅觉 相关 基因 数量 占 比 (%) 嗅觉 相关 基因 数量 占 比 (%) 
Number of olfaction-related genes Proportion Number of olfaction-related genes Proportion 

气味 结合 蛋白 基因 
US 1 3.70 5 13.16 
Odorant binding protein ( OBP) genes 
化 学 感受 蛋白 基因 

13 48. 15 17 44.74 
Chemosensory protein ( CSP) genes 
气味 / 响 觉 受 体 基因 

4 14.81 4 10.53 
Odorant/Olfactory receptor (OR) genes 
离子 型 受 体 基 因 

0 0 1 2.63 
Ionotropic receptor (IR) genes 
味觉 受 体 基因 

7 25.93 5 13.16 
Gustatory receptor ( GR) genes 
感觉 神经 元 膜 蛋白 基因 

2 7.41 6 15. 79 








Sensory neuron membrane protein (SNMP) genes 


MM 


3 it 
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迄今 为 止 ,嗅觉 相关 基因 已 经 在 20 RARA H 











( Theissinger et al., 2016) 比较 ,本 研究 获得 的 七 星 
IUE foh fA Fe oe fH P, se SE unigenes 占 比 为 
83. 2296 ,证 实 了 我 们 得 到 高 质量 的 七 星际 虫 触角 转 
录 组 测序 及 组 装 结果 。 采 用 NR 蛋白 数据 库 比 对 获 

















昆虫 中 被 研究 ,但 是 在 肚 虫 中 还 鲜 有 报道 。 此 外 , 通 
过 与 其 他 组 装 的 45 种 节肢 动物 转录 组 BUSCO 结果 


HAN 














得 的 unigenes 最 佳 比 对 物种 信息 注释 分 类 显示 ,最 
高 相似 性 的 物种 为 赤 拟 谷 资 和 中 欧 山 松 大 小 圳 (图 
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2) ,这 两 种 甲虫 作为 最 佳 比 对 物种 的 unigenes 占据 


2016; Hu et al., 2016b; Wang J et al., 2017) ; 3H 








了 所 有 注释 基因 的 约 70% ,这 主要 归 因 于 其 与 七 星 
球 虫 一 道 作为 鞘翅 目 昆 虫 , 具 有 相对 完整 的 全 基因 
组 序列 ,上 且 具 有 丰富 的 基因 组 注释 信息 ;尤其 赤 拟 谷 
盗 为 著名 的 模式 动物 ,拥有 鞘翅 目 昆 虫 中 最 全 面 的 
基 组 学 信息 (Tribolium Genome Sequencing 
Consortium, 2008 )。 已 有 大 量 研究 结果 表明 ,注释 
到 的 基因 均 与 基因 组 信息 已 知 且 亲缘 关系 较 近 的 物 
种 相似 性 序列 最 多 ( 郑海霞 等 , 2018; 杜 迎 刚 等 ， 
2019), GO 分 类 表明 ,七 星球 虫 触 角 转 录 组 基因 划 
分 为 生物 学 过 程 细胞 组 分 与 分 子 功能 3 个 亚 类 ,其 
中 参与 细胞 进程 .代谢 过 程 细胞 和 细胞 部 分 结合 、 
众 化 活性 的 基因 功能 注释 较 多 ,这 与 大 多 数 靖 翅 日 
昆虫 触角 转录 组 中 报道 的 (Cu et al., 2015; 郑海霞 
等 , 2018 ) 较为 一 致 ,表明 这 些 基 因 是 维持 甲虫 触角 
生物 学 功能 的 基础 。 当 然 , 不 同文 献 报道 的 不 同 鞘 



































目 美国 白 峨 .草地 量 Loxostege sticticalis ILE ARN 2k 
iR Eogystia hippophaecolus 触角 转录 组 中 分 别 为 124 , 
131 和 137 个 (Hu et al., 2016a; Zhang et al., 2016; 
Wei et al., 2017) ;2Ex8 ELS AN RC Diaphorina citri, 
TNR Subpsaltria yangi FU Œ JA £& Ej "f$. Cyrtorhinus 
lividipennis 触角 转录 组 中 分 别 为 126, 77 和 57 个 
(Wu et al., 2016; Qi et al., 2018; Wang et al., 
2018) ; 双 妈 上 日 南亚 实 刚 Bactrocera tau Wi A E he 
Episyrphus balteatus 4 K JK [b E $F "ig Eupeodes 
corollae WIAR I Bactrocera minax 触角 转录 组 中 
分 别 为 104, 154, 134 和 127 个 (Wang B et al., 
2017; 杜 迎 刚 等 , 2019; Wu et al., 2019) 。 因 此 ,我 
们 提出 昆虫 触角 嗅觉 相关 基因 的 绝对 数量 可 能 不 是 
昆虫 食性 分 化 一 个 最 重要 的 因素 。 如 两 种 食性 球 
E 不同 食 性 蝇 类 间 等 食性 差异 很 大 ,但 嗅觉 相关 基 




















翅 目 昆虫 触角 转录 组 中 基因 功能 注释 的 GO 数目 不 
[n] ,这 可 能 是 由 不 同 的 昆虫 种 类 、 生 存 环境 测序 手 
段 和 分 析 参 数 等 不 同 所 导致 的 。 此 外 ,包括 KEGG 
注释 信息 ,69.71% 的 unigenes 在 所 有 六 大 数据 库 中 
得 以 注释 ,这 个 结果 表明 有 超过 30% 的 七 星球 虫 触 
角 转 录 组 unigenes 没有 任何 注释 信息 , 可见 仍 有 一 
大 批示 知 基因 功能 的 新 基因 有 待 鉴定 ,其 中 也 可 能 
包含 许多 潜在 的 新 的 嗅觉 相关 基因 。 此 外 ， unigenes 
主要 涉及 代谢 、 免 疫 防 御 、 核 糖 体 和 和 蛋白 过 程 等 代谢 
通路 ,表明 了 外 源 物质 防御 、 免 疫 响 应 、 蛋 白质 的 合 
成 及 代谢 在 七 星 标 虫 嗅觉 行为 活动 中 的 重要 性 。 
本 研究 首次 以 两 种 食性 分 化 程度 最 高 的 球 虫 转 
录 组 为 对 象 , 从 七 星球 虫 中 挖掘 到 27 个 嗅觉 相关 基 
因 ( 表 2) ,并 通过 以 往 获 得 的 茄 二 十 八 星 球 虫 转录 
组 数据 (Zhang et al., 2018) ,从 中 鉴定 到 38 个 嗅觉 
相关 基因 ,主要 包括 OBP, CSP, OR, IR, GR 和 
SNMP 基因 。 此 外 ,已 有 研究 发 现 捕食 性 的 异 色 球 
Hi Harmonia axyridis 触角 转录 组 中 至 少 存在 57 个 
嗅觉 相关 基因 (包括 30 个 OBP 基因 和 27 个 CSP 基 












































因数 量 差异 较 小 ;相反 ,有 的 昆虫 类 群 (如 柑橘 木 虱 
与 枯 蝉 和 黑 肩 绿 盲 暑 、 红 脂 大 小 才 与 天 牛 和 象 甲 
等 ) 食 性 相似 ,但 嗅觉 相关 基因 数量 差异 较 大 。 另 
外 ,上 述 例子 中 两 种 ( 星 天 牛 和 光 肩 星 天 牛 ) 同 隶属 
Anoplophora 属 食 性 相似 的 天 牛 触角 转录 组 中 嗅觉 
相关 基因 数目 也 有 一 定 差异 (151 w 108 ) ,而 与 其 他 
畏 远 目 昆 虫 又 较为 一 致 ,这 两 种 天 牛 嗅 觉 相关 基因 
数量 差异 ,可 能 是 除 食性 外 的 其 他 环境 因子 引起 的 。 

在 七 星球 虫 甬 角 和 茄 二 十 八 星 球 虫 转录 组 中 ， 
CSP 基因 数量 均 占 据 了 嗅觉 相关 基因 数目 的 绝 大 多 
数 (40% 以 上 )( 表 2)。 不 同 昆虫 间 CSP 基因 所 占 
嗅觉 相关 基因 比重 不 同 ,如 星 天 牛 (10.60% ) 、 红 脂 
KU NIE (10% ) P H H IR (13.71% ) (Gu et al., 
2015; Zhang et al., 2016; Wang J et al., 2017) , 35 
外 ,OBP NEEESE f £f (3.70% ) 和 茄 二 十 八 
EUR (13. 16% ) 转录 组 中 所 占 比 例 差异 较 大 ( 表 
2) 。 其 他 嗅觉 相关 基因 类 型 也 是 如 此 ,可 见 不 同 昆 
虫 类 和 群 间 ,嗅觉 相关 基因 组 成 比重 不 尽 相 同 。 另 外 ， 
OBPs 是 嗅觉 发 生 进程 中 的 一 个 关键 环节 ,因为 它们 





























I) CK, 2018) 。 可 见 ,嗅觉 相 关 基 因数 量 差异 
在 捕食 -捕食 购 虫 间 比 捕食 - 植 食 又 虫 间 大 。 男 一 方 
面 ,不 同 食性 的 球 虫 间 ,嗅觉 相关 基因 数量 差异 不 及 
疆 虫 与 其 他 目 昆 虫 间 嗅觉 相关 基因 数量 差异 大 。 例 
JS H 21S Kha Dendroctonus valens 、 星 天 牛 
Anoplophora chinensis、 光 肩 星 天 和 牛 Anoplophora 
glabripennis 和 红 棕 象 甲 Rhynchophorus ferrugineus fih 
角 转 录 组 中 六 大 类 嗅觉 相关 基因 总 数 分 别 为 60， 
151, 108 和 157 个 (Gu et al., 2015; Antony et al., 





通过 昆虫 感觉 淋巴 将 气味 从 外 部 环境 传送 到 ORs 
(Leal, 2004, 2013) 。 在 七 星球 虫 触角 转录 组 中 ,未 
有 TIR 基因 被 鉴定 到 , 茄 二 十 八 星 飘 虫 中 鉴定 到 了 I 
个 IR 基因 ( 表 2)。IR 基因 家 族 包含 两 大 类 :一 类 
是 与 嗅觉 有 关 的 保守 的 触角 IRs; 男 一 类 是 在 包括 
味觉 器 官 在 内 的 其 他 组 织 中 表达 的 具有 物种 特异 性 
的 IRs ( Croset et al., 2010) 。 本 研究 中 ,七星 昧 虫 触 
角 转 录 组 未 发 现 IR 基因 ,再 者 茄 二 十 八 星 飘 虫 因为 
是 其 不 同 发 育 阶段 整个 虫 体 的 转录 组 混 池 组 装 结 
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果 , 故 这 些 证 据 表 明 球 虫 中 鉴定 到 的 TR 基因 极 有 可 
能 属于 第 2 种 类 型 。 然 而 ,还 需要 通过 与 其 他 昆虫 
两 种 类 型 的 IR. 基因 联合 构建 系统 发 育 树 和 进行 不 
同 组 织 基因 表达 分 析 最 终 确 认 。GR 基因 在 昆虫 中 
主要 在 味觉 带 官 的 味觉 感受 神经 元 中 表达 ,作为 糖 
受 体 、 热 传感器 二氧化碳 和 苦味 受 体 ,参与 糖 .苦味 
化 合 物 和 其 他 昆虫 信息 素 的 检测 (Ni et al., 2013; 
Gu et al., 2015) 。 本 研究 中 发 现 七 星 标 虫 触角 所 会 
GR 基因 的 比例 (25.93% ) 是 茄 二 十 八 星 标 虫 中 的 
W 2 倍 (13.16% )( 表 2) ,表明 了 七 星 味 虫 在 捕食 过 
程 中 触角 有 更 强 的 代谢 能 力 , 对 糖 及 昔 味 化 合 物 有 
更 灵敏 的 嗅觉 ,对 其 他 昆虫 (尤其 捕食 对 象 ) 有 更 快 
速 的 响应 。SNMPs 的 功能 在 昆虫 中 还 鲜 有 报道 , 故 
仍 需 更 多 的 实验 研究 。 

本 研究 测序 .组 装 及 注释 了 七 星 标 虫 触角 转录 
组 ,并 结合 茄 二 十 八 星 标 虫 不 同 发 育 阶段 测序 数据 
混 池 组 装 的 转录 组 ,在 两 种 食性 的 又 虫 中 鉴定 了 六 
大 类 昆虫 嗅觉 相关 基因 ,并 基于 每 类 基因 所 占 鉴定 
出 的 嗅觉 相关 基因 总 数 比 例 ,比较 了 它们 在 两 种 食 
性 又 虫 间 的 差异 。 值 得 注意 的 是 ,因为 测序 深度 及 
测序 组 织 的 不 同 ,直接 比较 每 种 嗅觉 相关 基因 类 型 
在 两 种 飘 虫 中 的 绝对 数目 是 从 忆 的 ,故我 们 将 其 转 
换 成 统一 测序 背景 下 ,不 同类 型 嗅觉 相关 基因 在 总 
数 中 所 占 的 比重 /比例 进行 比较 。 本 研究 获得 了 七 
星 标 虫 高 质量 的 触角 转录 组 及 两 种 球 虫 中 的 一 批 嗅 
觉 相关 基因 ,但 这 些 基因 的 组 织 表达 模式 和 生理 生 
态 功能 仍 有 待 进一步 研究 。 今 后 ,可 利用 实时 荧光 
定量 POR 技术 ,对 早 虫 嗅觉 相关 基因 组 织 特异 表达 
情况 确认 ,以 辅助 明确 其 生物 学 功能 ;结合 不 同 层次 
的 组 学 测序 结果 ,利用 正 反 向 遗传 学 手段 ,过 表达 或 
沉默 关键 嗅觉 相关 基因 ,进一步 明确 嗅觉 相关 基 
TESI dU E TET EP WE HI s 
































参考 文献 (References) 


Antony B, Soffan A, Jakše J, Abdelazim MM, Aldosari SA, Aldawood 
AS, Pain A, 2016. Identification of the genes involved in odorant 
reception and detection in the palm weevil Rhynchophorus 
ferrugineus , an important quarantine pest, by antennal transcriptome 
analysis. BMC Genomics, 17 ; 69. 

Bolger AM, Lohse M, Usadel B, 2014. Trimmomatic; a flexible trimmer 
for Illumina sequence data. Bioinformatics, 30(15) : 2114 —2120. 

Brown EB, Layne JE, Zhu C, Jegga AG, Rollmann SM, 2013. Genome- 
wide association mapping of natural variation in odor-guided behavior 
in Drosophila. Genes Brain Behav., 12(5) : 503 —515. 

Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, 


Gibson TJ, Benton R, Stern DL, 2010. Ancient protostome origin 
of chemosensory lonotropic glutamate receptors and the evolution of 
insect taste and olfaction. PLoS Genet., 6(8) : e1001064. 

Du YG, Ji QE, Lai ZX, 2019. Preliminary analysis of transcriptome and 
olfaction-related genes in Bactrocera tau ( Walker). Chin. J. Trop. 
Crops, 40(9) : 1796 - 1803. [ 杜 迎 刚 , 季 清 娥 ， 赖 钟 雄 ，2019. 
南亚 实 蝇 转 录 组 及 嗅觉 相关 基因 的 初步 分 析 . 热带 作物 学 报 ， 
40(9) : 1796 -1803 ] 

Escalona HE, Zwick A, Li HS, Li J, Wang X, Pang H, Hartley D, 
Jermiin LS, Nedvěd O, Misof B, Niehuis O, Slipifiski A, 
Tomaszewska W, 2017. Molecular phylogeny reveals food plasticity 
in the evolution of true ladybird beetles ( Coleoptera; Coccinellidae : 
Coccinellini). BMC Evol. Biol., 17(1) : 151. 

Field LM, Pickett JA, Wadhams LJ, 2001. Molecular studies in insect 
olfaction. Insect Mol. Biol., 9(6) : 545 —551. 

Gadenne C, Barrozo RB, Anton S, 2016. Plasticity in insect olfaction : 
to smell or not to smell? Annu. Rev. Entomol., 61: 317 —333. 

Gu XC, Zhang YN, Kang K, Dong SL, Zhang LW, 2015. Antennal 
transcriptome analysis of odorant reception. genes in the red 
turpentine beetle ( RTB ), Dendroctonus valens. PLoS ONE, 10 
(5) : e0125159. 

Hu P, Tao J, Cui M, Gao C, Lu P, Luo Y, 2016a. Antennal 
transcriptome analysis and expression profiles of odorant binding 
proteins in Eogystia hippophaecolus ( Lepidoptera: Cossidae). BMC 
Genomics, 17 ; 651. 

Hu P, Wang J, Cui M, Tao J, Luo Y, 2016b. Antennal transcriptome 
analysis of the Asian longhorned beetle Anoplophora glabripennis. 
Sci. Rep., 6; 26652. 

Jia XJ, Wang HX, Yan ZG, Zhang MZ, Wei CH, Qin XC, Ji WR, 
Falabella P, Du YL, 2016. Antennal transcriptome and differential 
expression of olfactory genes in the yellow peach moth, Conogethes 
punctiferalis ( Lepidoptera; Crambidae). Sci. Rep., 6; 29067. 

Leal WS, 2004. Pheromone reception. Top. Curr. Chem., 240; 341 — 
360. 

Leal WS, 2013. Odorant reception in insects; roles of receptors, binding 
proteins, and degrading enzymes. Annu. Rev. Entomol., 58; 373 — 
391. 

Matsubayashi KW, Kohyama TI, Kobayashi N, Yamasaki S, Kuwajima 
M, Katakura H, 2017. Genetic divergence with ongoing gene flow is 
maintained by the use of different hosts in phytophagous ladybird 
beetles genus Henosepilachna. J. Evol. Biol., 30 (6): 1110 — 
1123. 

Ni L, Bronk P, Chang EC, Lowell AM, Flam JO, Panzano VC, 
Theobald DL, Griffith LC, Garrity PA, 2013. A gustatory receptor 
paralogue controls rapid warmth avoidance in Drosophila. Nature, 
500(7464) : 580 —584. 

Pang H, 2002. Current situation of the systematics of Coccinellidae. 
Entomol. Knowl., 39(1): 17 -22. [ 庞 虹 , 2002. URPI E 
究 的 现状 . 昆虫 知识 , 39(1) ; 17-22] 

Qi M, Wei S, Wei C, 2018. Identification of candidate olfactory genes in 








cicada Subpsaltria yangi by antennal transcriptome analysis. Comp. 


Biochem. Physiol. D Genomics Proteomics, 28; 122 —133. 


726 昆虫 学 报 Acta Entomologica. Sinica 63 35 





Saito K, Nomura S, Yamamoto S, Niiyama R, Okabe Y, 2017. 
Investigation of hindwing folding in ladybird beetles by artificial 
elytron transplantation and microcomputed tomography. Proc. Natl. 
Acad. Sci. USA, 114(22) : 5624 —5628. 

Sato K, Touhara K, 2009. Insect olfaction; receptors, signal 
transduction, and behavior. Results Probl. Cell Differ., 47; 121 — 
138. 

Tang Q, Zang G, Cheng C, Luan M, Dai Z, Ying X, Yang Z, Zhao L, 
Su J, 2017. Diplosporous development in Boehmeria tricuspis ; 
insights from de novo transcriptome assembly and comprehensive 
expression profiling. Sci. Rep., 7; 46043. 

Theissinger K, Falckenhayn C, Blande D, Toljamo A, Gutekunst J, 
Makkonen J, Jussila J, Lyko F, Schrimpf A, Schulz R, Kokkoc H, 
2016. De Novo assembly and annotation of the freshwater crayfish 
Astacus astacus transcriptome. Mar. Genomics, 28; 7 —10. 

Tribolium Genome Sequencing Consortium, 2008. The genome of the 
model beetle and pest Tribolium castaneum. Nature, 452 (7190) : 
949 - 955. 

Van der Werf W, Evans EW, Powell JA, 2000. Measuring and 
modelling the dispersal of Coccinella septempunctata ( Coleoptera: 
Coccinellidae) in alfalfa fields. Eur. J. Entomol., 97 (4): 487 — 
493. 

Wang B, Liu Y, Wang GR, 2017. Chemosensory genes in the antennal 
transcriptome of two syrphid species, Episyrphus balteatus and 
Eupeodes corollae ( Diptera: Syrphidae). BMC Genomics, 18 (1): 
586. 

Wang GY, Zhu JL, Zhou WW, Liu S, Khairul QM, Ansari NA, Zhu 
ZR, 2018. Identification and expression analysis of putative 
chemoreception genes from Cyrtorhinus lividipennis ( Hemiptera: 
Miridae) antennal transcriptome. Sci. Rep., 8: 12981. 

Wang J, Hu P, Gao P, Tao J, Luo Y, 2017. Antennal transcriptome 
analysis and expression profiles of olfactory genes in Anoplophora 
chinensis. Sci. Rep., 7; 15470. 

Wei HS, Li KB, Zhang S, Cao YZ, Yin J, 2017. Identification of 
candidate chemosensory genes by transcriptome analysis in Loxostege 


sticticalis Linnaeus. PLoS ONE, 12(4) : e0174036. 


Wu Z, Kang C, Qu M, Chen J, Chen M, Bin S, Lin J, 2019. 
Candidates for chemosensory genes identified in the Chinese citrus 
fly, Bactrocera minax, through a transcriptomic analysis. BMC 
Genomics , 20(1) : 646. 

Wu Z, Zhang H, Bin S, Chen L, Han Q, Lin J, 2016. Antennal and 
abdominal transcriptomes reveal chemosensory genes in the Asian 
citrus psyllid, Diaphorina citri. PLoS ONE, 11(7) : e0159372. 

Zhang LW, Kang K, Jiang SC, Zhang YN, Wang TT, Zhang J, Sun L, 
Yang YQ, Huang CC, Jiang LY, Ding DG, 2016. Analysis of the 
antennal transcriptome and insights into olfactory genes in 
Hyphantria cunea ( Drury). PLoS ONE, 11; e0164729. 

Zhang QL, Wang F, Guo J, Deng XY, Chen JY, Lin LB, 2018. 
Characterization of ladybird Henosepilachna vigintioctopunctata 
transcriptomes across various life stages. Sci. Data, 5; 180093. 

Zhang S, Pang B, Zhang L, 2015. Novel odorant-binding proteins and 
their expression patterns in grasshopper,  Oedaleus asiaticus. 
Biochem. Biophys. Res. Commun., 460(2) : 274 —280. 

Zhang S, Zhang Z, Wang H, Kong X, 2014. Antennal transcriptome 
analysis and comparison of olfactory genes in two sympatric 
defoliators, ^ Dendrolimus houi and Dendrolimus — kikuchii 
(Lepidoptera: Lasiocampidae ). Insect Biochem. Mol. Biol., 52 
(1): 69-81. 

Zhang W, 2018. Identification of Preference Volatiles to Harmonia 
axyridis in the Aphid-infested Plants and Exploration of Olfactory 
Related Genes. MSc Thesis, Yangzhou University, Yangzhou, 
Jiangsu. [ 张 伟 , 2018. $ k fé E REL] aer A AUR P FER W B6] 
鉴定 及 其 嗅觉 相关 基因 的 挖掘 . 江苏 扬州 : 扬州 大 学 硕士 学 位 
论文 ] 

Zheng HX, Zhang YW, Zhang XH, Yang MH, 2018. Analysis of the 


























antennal transcriptome and olfaction-related genes of Callosobruchus 
chinensis ( Coleoptera; Bruchuidae). Acta Entomol. Sin., 61(2): 
168 -177. [郑海霞 , KE oc, 张 仙 红 , 杨 美 红 , 2018. 绿豆 象 
触角 转录 组 及 嗅觉 相关 基因 的 分 析 . 昆虫 学 报 , 61(2) : 168 - 
177] 























(责任 编辑 : 马 丽 萍 ) 


